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escence, and  IRFET  measurements  were  performed  to  characterize  the 
X level  and  to  attempt  to  determine  its  nature. 
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SECTION  1 


INTRODUCTION 

The  objective  of  the  silicon  indium  material  and  detector 
characterization  (SIMDEC)  program  being  conducted  by  Hughes  under  con- 
tract DAAK70-77-C-0082  is  to  investigate  the  nature  of  defects  in  indium- 
doped  silicon  detector  material  and  to  develop  techniques  for  the  control 
of  quality  in  the  growth  of  extrinsic  and  high-purity  undoped  silicon  for 
infrared  (IR)  detector  monolitic  chips.  The  discovery  of  a second,  shal- 
lower acceptor  level  in  Si: In,  the  so-called  X level,  demonstrated  the 
need  for  further  research  on  indium-doped  silicon.  The  X level  has  been 
observed  in  measurements  of  Hall  effect  versus  temperature  and  of  IR 
photoconductivity  and  absorption.  It  is  undesirable  in  extrinsic  Si:In 
photoconductors  because  its  shallower  energy  level  results  in  excess 
thermally  ionized  carriers.  These  lower  the  maximum  temperature  at  which 
background-limited  detectivities  can  be  obtained.  This  can  reduce  the 
operating  temperature  by  as  much  as  10  K with  a corresponding  increase 
in  cooling  power  requirements. 

The  X level  has  an  ionization  energy  of  0.11  eV  and  has  a concen- 
tration that  is  a direct  function  of  the  concentration  of  the  indium. 

Different  growth  conditions  can  produce  drastically  different  amounts  of 
X level  in  comparably  indium  doped  material  samples.  Thus,  it  cannot  be 
a discrete  dopant.  These  facts  indicate  that  the  X level  is  likely  to 

, 

be  a complex  of  indium  and  some  other  crystalline  defect  or  impurity. 

However,  the  exact  physical  nature  of  the  X level  is  unknown,  as  are 
the  growth  and  processing  factors  that  determine  its  concentration.  The 
purpose,  then,  of  the  SIMDEC  program  is  to  apply  the  appropriate  physi- 
cal measurements  to  Sirin  to  identify  the  nature  of  the  X level  and  means 
of  controlling  it. 

The  prime  contractor  for  this  program  is  Hughes  Research  Laboratories 
(HRL);  the  California  Institute  of  Technology  (CIT)  and  the  University  of 
California  it  Davis  (UCD)  are  the  subcontractors. 


PRECEDING  PAGE  BLANK. NOT  FILMED 


The  most  significant  result  during  this  reporting  period  is  the 
discovery  of  a correlation  between  the  X-level  concentration  and  the 
concentration  of  carbon.  This  relationship  was  determined  from  X-level 
data  (taken  by  Hall-effect  measurements)  and  C data  (from  IR  transmission 
studies) . This  C dependence  explains  the  higher  than  expected  X-level 
concentrations  observed  in  float-zone  samples,  which  were  found  to  have 
C concentrations  higher  than  expected.  It  also  explains  the  lower  X- 
level  values  in  samples  from  the  large  Czochralski  grower,  which  were 
found  to  have  low  C. 

Although  these  results  must  be  considered  preliminary  pending  data 
on  more  crystals,  a model  for  the  X level  can  be  proposed.  The  X level 
may  be  an  In-C  pair  occupying  adjacent  substitutional  sites  in  the 
lattice.  The  lower  X-level  concentrations  observed  after  850°C  anneals 
(as  compared  with  650°C  anneals)  are  consistent  with  this  model  if  a 
reasonable  In-C  pair  binding  energy  is  assumed.  There  will  be  fewer 
In-C  pairs  in  thermal  equilibrium  at  850°C  than  at  650°C  because  of  the 
higher  lattice  thermal  energy.  This  concentration  of  pairs  is  "frozen" 
by  the  quench  following  anneal.  A subsequent  anneal  at  650°C  will 
increase  the  number  of  pairs,  which  is  consistent  with  the  observed 
higher  X-level  concentration  at  this  anneal  temperature.  It  is  not 
necessary  to  invoke  Yoshihiro’ s'*'  observation  of  a minimum  in  substitu- 
tional carbon  at  800°C,  which  in  fact  we  do  not  see  with  our  much 
shorter  anneals. 

The  source  of  the  carbon  in  the  Czochralski-grown  Si: In  is  probably 
the  graphite  heaters  used  to  melt  the  poly-silicon.  Since  the  heater  is 
not  used  at  full  power  in  the  large  Czochralski  grower,  less  C is  likely 
to  be  vaporized  from  it  during  growth.  Since  a causal  relationship 
between  C and  the  X level  has  been  found,  it  now  becomes  possible  to 
consider  ways  to  modify  the  Czochralski  growers  to  minimize  the  X-level 
concentration  and  eliminate  this  problem  for  Sirin  photoconductive 
detectors . 
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The  results  of  studies  carried  out  at  HRL  are  contained  in  Section  2. 
Two  additional  crystals  were  grown  for  the  program  during  this  period, 
and  the  results  of  Hall-effect  analysis  are  presented  and  compared  with 
the  predicted  dopant  profiles.  IR  transmission  data  and  the  concentra- 
tions of  indium,  oxygen,  and  carbon  obtained  are  discussed.  Also  discussed 
are  the  Hall-effect  evaluations  of  material  processed  at  UCD. 

An  intentionally  uncompensated  crystal  was  found  to  contain  an 
excessive  amount  of  boron;  the  origin  of  the  boron  is  presently  being 
investigated.  The  second  crystal  was  a simulation  of  Si: In  crystals 
pulled  at  Hughes  Industrial  Products  Division,  Carlsbad,  California  (IPD). 
Since  this  crystal  contained  a high  X-level  concentration,  several  growth 
parameters  (e.g.,  argon  overpressure,  melt  surface  area  to  volume  ratio, 
crucible  material,  and  silicon  starting  material)  were  eliminated  as 
potential  causes  of  the  X level.  The  spectroscopic  studies  indicated 
that  the  X level  is  uncorrelated  to  the  oxygen  concentration  but  is 
strongly  correlated  to  the  carbon  concentration  in  the  crystal.  Several 
changes  were  observed  in  the  material  processed  at  UCD,  indicating  that 
dopant  concentrations  obtained  on  the  original  ingot  should  not  be  relied 
on  when  making  device  calculations. 

The  progress  at  CIT  during  this  period  is  described  in  Section  3. 

The  three  strongest  photoluminescence  lines  observed  uniquely  in  Sirin 
were  investigated.  The  temperature  and  time  dependencies  of  these  line 
strengths  were  measured.  The  results  indicate  that  these  lines  are  due 
to  an  isoelectronic  trap  (or  traps)  and  not  to  an  acceptor,  such  as  the 
X level. 

Characterizations  of  deep  level  impurity  center  concentrations, 
energy  levels,  thermal  and  optical  emission  rates,  and  thermal  capture 
cross  sections  were  performed  by  UCD  during  this  reporting  period  using 
the  MOSFET  structure.  Indium  concentrations  measured  by  several  methods 
(variation  of  the  MOSFET  threshold  voltage  with  substrate  bias,  reverse 
bias  junction  capacitance,  and  resistivity  versus  temperature)  were  found 
to  be  typically  a factor  of  two  less  than  results  from  Hall-effect  mea- 
surements. The  indium  hole  emission  rate  at  77  K was  found  to  be  6.0  msec 
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with  a 0.15  cV  activation  energy  for  a lightly  doped  sample  and  20  ysec 

and  0.117  eV  for  a more  heavily  doped  sample.  The  hole  capture  coefficient 

-7  3 

was  found  to  be  3.7  x 10  cm  /sec.  Small  transients  in  the  thermal 

emission  rate,  associated  with  the  X level,  were  observed,  which  indi- 

“6  3 

cate  a capture  coefficient  of  1.5  x 10  cm  /sec  at  35  K for  this  level. 
Measurements  of  Si: In  photoconductors  are  also  included.  The  description 
of  this  work  will  be  published  later. 

The  appendices  to  this  report  contain  reprints  of  papers  given  at 
the  1978  IRIS  Detector  Specialties  Group  Meeting. 
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SECTION  2 


CRYSTAL  GROWTH  AND  CHARACTERIZATION 


This  section  describes  progress  at  HRL  during  this  reporting  period. 
The  growth  and  analysis  of  two  additional  crystals  are  discussed.  An 
evaluation  program  to  characterize  silicon  in  terms  of  its  oxygen  and 
carbon  content  using  IR  spectroscopy  is  described,  and  the  initial 
results  are  presented.  Results  of  Hall-effect  analysis  of  material 
processed  at  UCD  are  also  given. 

A.  CRYSTAL  GROWTH 


Ingot  C076  was  grown  in  a 1-kg  system  as  part  of  the  SIMDEC  program 
to  produce  a very  lightly  counterdoped  crystal  with  a low  boron  concen- 
tration. The  principal  sources  of  boron  in  Czochralski-grown  silicon  are 
the  starting  polycrystalline  silicon  and  the  quartz  crucible,  which  dis- 
solves slowly  in  the  melt  during  growth  and  releases  its  boron  impurity. 
Low  boron  starting  material  was  used  and  the  growth  was  from  a Suprasil- 
lined  crucible.  This  type  of  container  was  used  because  Suprasil  has  a 
significantly  lower  boron  content  than  does  regular  quartz,  which  is 
typically  used  for  crucibles.  Since  B cannot  be  removed  from  Si  during 
the  growth  process  to  any  significant  extent,  care  must  be  taken  to  mini- 
mize its  introduction  during  growth.  Although  Suprasil  has  a higher  dis- 
solution rate  into  Si,  the  B buildup  in  the  melt  is  considerably  slower 
than  for  regular  quartz.  This  results  in  a crystal  with  lower  B 
concentration. 

Since  P must  be  added  to  compensate  for  the  B,  crystals  grown  with 
this  type  of  crucible  require  less  P to  be  added  to  ensure  that  Np  > Ng 
throughout  the  ingot.  Therefore,  they  will  have  a lower  net  compensation 
and  a corresponding  higher  detector  responsivity . Since  the  values  of 
Np  and  Ng  are  lower,  their  difference  can  also  be  kept  smaller. 

No  P was  intentionally  added  during  the  growth  of  crystal  C076. 

The  major  sources  of  P were  the  silicon  starting  material,  which  had 

12  -i 

an  initial  P concentration  of  8 x 10  cm" J,  and  the  crucible.  Phosphorus 


tends  to  build  up  in  the  melt  (and  the  crystal)  during  growth  because 
its  segregation  coefficient  is  less  than  one.  The  rate  of  this  buildup 
is  reduced  by  evaporation  of  P from  the  melt  under  the  low-pressure  con- 
dition of  the  crystal  growth.  The  overall  result  under  these  conditions 
is  a net  increase  in  B relative  to  P in  the  crystal  from  seed  to  tang  due 
to  the  rate  of  B buildup  from  crucible  dissolution. 

Table  1 lists  the  results  of  Hall-effect  versus  temperature 
measurements  on  three  samples  from  ingot  C076.  The  samples  labeled  "A" 
are  from  the  seed  end  of  the  crystal.  The  "M"  sample  is  from  about  mid- 
way between  the  seed  and  tang.  All  samples  are  from  the  center  of  the 
wafer  from  which  they  were  cut.  Sample  A20  received  a 650°C  anneal  for 
30  min  to  reduce  the  oxygen  donor  content.  Samples  A30  and  M30  received 
an  850°C  anneal  for  60  min,  which  we  have  found  reduces  the  X-level 
concentration  and  the  oxygen-donor  concentration. 

The  850°C  anneal  reduced  the  X level  as  can  be  seen  by  comparing 
samples  A20  and  A30.  In  addition,  the  X level  increases  in  the  crystal 
from  seed  to  tang.  Samples  A20  and  A30,  which  were  adjacent,  differ  by 
an  order  of  magnitude  in  net  compensation,  N^-N^.  This  is  apparently 
due  to  slight  fluctuations  in  an  extremely  close  balance  between  the 
P and  B concentrations.  Both  B and  P increased  in  concentration  toward 
the  center,  axially,  of  the  ingot;  B increased  more  rapidly,  as  expected. 


Table  1.  Dopant  Concentrations  in  Crystal  C076 


An  analysis  of  the  pickup  of  B and  P in  crystal  C076  was  made  using 

2 

a model  developed  by  H.  Kimura.  The  dopant  concentrations  measured  on 

the  Hall  samples  were  used  as  boundary  conditions,  and  initial  values  of 

14  3 14 

B and  P in  the  melt  were  calculated  to  be  1.7  x 10  B/cm  and  8 x 10 

3 

P/cm  . This  amount  of  B far  exceeds  the  amount  possible  due  to  crucible 
dissolution.  However,  it  corresponds  to  0.6  ppm  in  the  indium  material 
used  for  doping.  This  is  within  the  manufacturer's  specification  for 
purity  of  6-9 's  indium.  Other  Si: In  crystals  doped  with  this  material 
were  also  doped  with  P and  were  grown  in  regular  quartz  crucibles;  as  a 
result,  this  amount  of  B pickup  would  not  be  observed. 

The  results  of  measurements  of  crystal  C076  suggested  several 
experiments,  which  are  now  underway.  A section  of  this  crystal  is  under- 
going several  passes  of  vacuum  float  zoning  to  remove  the  In  and  P, 
leaving  only  the  B in  its  original  distribution.  This  will  allow  a 
verification  of  the  model  on  which  the  calculations  were  based.  Another 
crystal  is  being  grown  using  In  obtained  from  IPD.  This  material  is 
known  to  have  a lower  B concentration.  The  In  dopant  material  currently 
being  used  will  be  baked  under  a controlled  atmosphere  to  reduce  its 
B content.  Samples  of  In  from  HRL  and  IPD  will  be  analyzed  for  impuri- 
ties by  mass  spectroscopic  techniques. 

Crystal  C055  was  an  experiment  to  simulate,  using  a small  (1  kg) 
grower  at  HRL,  growth  of  Si:In  crystals  in  the  large  crystal  grower  at 
IPD.  The  same  polysilicon  starting  material  and  In  and  P dopant  materials 
were  used,  and  the  gas  overpressure  and  the  pull  and  rotation  rates  were 
the  same.  The  ratio  of  crystal  diameter  to  crucible  diameter  was  also 
identical.  The  object  of  this  simulation  was  to  determine  if  any  of 
these  growth  parameters  influenced  the  X level.  The  results  of  Hall- 
effect  measurements  of  samples  from  the  center  of  the  seed  end  of  this 
ingot  are  presented  in  Table  2.  Sample  A30  had  a 650°C  anneal  for 
30  min,  and  sample  A40  had  an  850°C  anneal  for  1 hr.  The  ingot  has  a 
rather  high  X-level  concentration,  one  similar  to  the  levels  observed 
in  crystals  grown  under  the  normal  conditions  for  this  machine.  Since 
this  crystal  was  grown  under  the  same  conditions  as  crystals  at  IPD  that 
consistently  exhibit  X-level  concentrations  at  or  below  10^  cm  \ it 
is  unlikely  that  these  growth  parameters  are  influencing  the  X level. 
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Table  2.  Hall-Effect  Results  for  the  Seed  End  of  Crystal  C055 


6303 


B.  IR  SPECTROSCOPIC  STUDIES 

IR  measurements  of  Si: In  grown  under  different  conditions  and  with 
varying  In  doping  levels  were  begun  during  this  reporting  period  to 
determine  the  oxygen  and  carbon  content  of  the  crystals  and  the  relation- 
ship, if  any,  of  these  impurities  to  the  X-level  concentration.  A 

correlation  between  In  concentration  and  X-level  concentration  has  been 

3*~5  6—8 

reported  and  substantiated  by  others.  However,  since  samples 

grown  under  different  conditions  with  similar  In  concentrations  show  dif- 
ferences in  the  X-level  concentration,  it  is  unlikely  that  the  X level 
is  a simple  In  complex.  At  least  one  more  component  appears  to  be 
necessary . 

The  experimental  apparatus  consisted  of  two  liquid-nitrogen  dewars 
designed  and  fabricated  at  HRL;  the  dewars  are  compact  enough  to  fit 
into  the  sample  and  reference  compartments  of  the  Beckman  IR-12  IR 
spectrometer.  The  double-beam-difference  method  was  used  to  obtain  the 
absorption  spectrum  of  the  sample.  An  oxygen-  and  carbon-free  sample  of 
undoped,  zone-refined  Si  (furnished  by  John  Baker,  Dow  Corning)  was 
mounted  in  the  reference  dewar,  and  the  sample  to  be  tested  was  mounted 
in  the  matching  dewar  in  the  sample  beam.  Both  dewars  have  Csl  windows 
that  allow  measurements  to  be  made  to  30  pm.  The  double-beam-difference 
method  canceled  the  lattice  absorptions  that  appear  near  both  the  oxygen 
and  carbon  impurity  center  peaks.  Because  of  the  lattice  absorption 
interference,  care  was  taken  to  assure  that  both  the  unknown  and  the 
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reference  samples  were  of  the  same  thickness.  (This  measurement 

technique  and  others  are  described  in  Refs.  9 and  10.)  There  has 

traditionally  been  disagreement  among  experimenters  as  to  the  correct 

calibration  factor  to  use  to  convert  optical  data  to  oxygen  and  carbon 

concentrations.^  The  relationships  used  in  this  study  between  the  peak 

absorption  coefficients  at  liquid  nitrogen  temperature  (an  and  a ) and 

3 u c 

concentrations  (atoms/cm  ) are 

[0]  = 1.43  x 1017  aQ 
[C]  = 3.1  x 1016  ac 

The  oxygen  calibration  is  derived  from  Refs.  9 and  10  scaled  to  77  K, 
and  the  carbon  calibration  is  from  ASTM  F123-74  corrected  for  a differ- 
ent baseline.  The  limits  of  sensitivity  for  measuring  oxygen  and  carbon 

15  3 15  3 

concentrations  are  1 x 10  atoms/cm  and  2 x 10  atoms/cm  , 

respectively. 

The  results  obtained  to  date  are  given  in  Table  3.  Additional 
optical  and  Hall-effect  measurements  are  in  progress  on  these  and  other 
samples.  Figure  1 shows  X- level  concentration  normalized  by  the  In  con- 
centration (both  obtained  from  Hall  measurements)  plotted  against  the 
interstitial  oxygen  concentration  obtained  from  the  optical  measurements. 

The  circles  are  HRL  measurements  and  are  labeled  with  the  crystal  number. 

■ 

The  triangles  are  data  from  solution-grown  Si: In  crystals  grown  and 

12 

measured  at  Honeywell.  Clearly,  no  correlation  exists  with  the  oxygen 
concentration.  Sample  Z008  is  particularly  significant  in  this  regard  as 
it  has  at  least  three  orders  of  magnitude  less  oxygen  and  yet  is  compara- 
ble to  the  other  samples  in  X level. 

In  contrast  to  the  lack  of  correlation  seen  for  the  X level  and 
oxygen  concentrations.  Figures  2 and  3 show  a definite  correlation  to 
carbon  concentration.  In  Figure  2,  the  X-level  concentration  normalized 
to  the  In  concentration  is  plotted  versus  the  carbon  impurity  concentra- 
tion. In  Figure  3,  the  X-level  concentration  is  plotted  against  the 
carbon  concentration. 

I 
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Summary  of  Hall-Effect  and  IR  Transmission  Measurements 
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C.  ANALYSIS  OF  UCD-PROCESSED  MATERIAL 

HRL,  in  support  of  its  subcontract  with  UCD,  performed  Hall-effect 
measurements  and  analyses  on  several  samples  that  had  been  processed  at 
UOD.  This  was  done  to  correlate  measurements  at  the  two  facilities  and 
to  provide  UCD  with  the  parameters  of  their  as-processed  samples  in  addi- 
tion to  the  as-grown  values. 

During  the  previous  period,  a Hall  sample  from  the  edge  of  a wafer 
processed  at  UCD  was  fabricated  and  measured.  Also,  a "bulk  resistor" 
manufactured  at  UCD  was  measured  at  HRL.  The  resistor  measurements 
agreed  with  those  at  UCD,  but  not  with  resistivity  data  taken  on  the 
Hall  sample  from  the  same  processed  wafer.  During  the  present  period, 
the  actual  device  was  fabricated  into  a Hall  sample  and  measured.  The 
results  are  shown  in  Figure  4,  which  is  a plot  of  resistance  versus 
temperature.  The  solid  line  and  discrete  points  are  the  bulk  resistor 
measurements  of  HRL  and  UCD,  respectively.  The  dashed  line  is  from  the 
edge  of  the  UCD-processed  wafer,  and  the  dash-dotted  line  is  from  the  / 
Hall  sample  made  from  the  bulk  resistor.  The  last  two  curves  were 


The  vertical  lines  for  samples  Z008  and  C076  represent  the 
uncertainty  in  X-level  concentration  for  these  samples.  This  uncertainty 
occurs  because  the  optical  and  Hall-effect  crystals  were  not  from  the 
same  portion  of  the  ingot.  There  is  also  uncertainty  concerning  the 
relative  calibrations  between  the  HRL  and  Honeywell  data.  The  data  were 
taken  at  different  temperatures  and  only  absorption  coefficients  were 
reported  by  Honeywell.  Thus,  while  relationships  to  the  carbon  concen- 
tration can  be  ascertained , relative  concentrations  cannot  be  obtained. 

The  carbon  concentrations  for  the  Honeywell  samples  represent  our  best 
estimate  of  the  appropriate  calibration  factor. 

It  appears  that  the  ratio  Nx/Njn  gives  the  better  correlation  to  the 
carbon  concentration,  but  there  are  insufficient  data  for  a conclusive 
result  to  be  obtained.  Other  samples  currently  are  being  measured;  these 
samples  represent  a wider  range  of  the  parameters. 
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obtained  by  multiplying  the  measured  resistivity  data  by  18,  which  is  the 
ratio  of  resistance  to  resistivity  of  these  devices  calculated  by  UCD. 

The  bulk  resistor  gives  a good  measure  of  the  substrate  resistivity, 
except  at  the  higher  temperatures.  This  error  is  apparently  due  to  contact 
resistance  in  the  bulk  resistor. 

Figure  4 indicates  a significant  difference  between  the  center  and 
edge  samples  from  this  wafer.  This  effect  is  also  seen  in  the  dopant 
data  from  Hall-effect  measurements,  which  is  given  in  Table  4.  Sample 
C010LF.10A  is  the  Hall  sample  from  the  bulk  resistor,  and  sample 
C010LF.410  is  from  the  edge  of  the  same  processed  wafer  in  a diffusion- 
free  region.  The  indium  concentration  is  not  significantly  different  in 
the  two  samples,  but  the  X level  is  dramatically  lower  and  the  net  donor 
concentration  is  also  lower.  This  result  is  not  consistent  with  data 
from  unprocessed  wafers,  which  do  not  show  large  variations  from  center 
to  edge. 

Table  4 also  contains  data  from  the  edge  of  another  UCD-processed 
wafer,  D,  from  crystal  C010  and  the  data  from  seed  and  tang  of  the 
original  ingot.  This  sample  showed  an  increase  in  net  donor  concentra- 
tion after  processing.  Data  from  the  edge  of  a UCD-processed  wafer  from 
crystal  C042  and  the  original  seed  and  tang  data  are  included  in  the 
table.  The  processed  edge  sample  again  shows  an  increase  in  net  donor 
concentration.  In  this  case,  the  X level  also  increased  slightly.  The 
specific  cause  of  these  processing  effects  is  unknown,  and  their  investi- 
gation is  beyond  the  scope  of  this  program. 
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Table  4.  Dopant  Concentrations  of  UCD-Processed 

Material 


Sample 

nin’ 

1015  cm-3 

NX’ 

1014  cm"3 

W 

1014  cm-3 

C010LF.D10 

Run  13  Edge 

1.3 

3.4 

3.7 

C010LF.10A 

Run  4 No.  10 

2.0 

0.37 

0.15 

C010LF.410 

Run  4 Edge 

2.2 

10 

0.47 

C010 

Seed-Tang 

0.95-4.0 

7.6-29 

0.86-0.44 

CO 42LF.10A 

Proc.  Edge 

2.3 

1.0 

1.4 

C042 

Seed-Tang 

1.9-2. 3 

0.29-0.32 

0.67-0.67 

6303 
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SECTION  3 

PHOTOLUMINESCENCE  STUDIES  IN  EXTRINSIC  SILICON 

This  section,  which  describes  the  work  done  at  CIT  during  this 
reporting  period,  was  written  by  Prof.  Thomas  McGill. 

The  primary  aim  of  the  work  has  been  to  explore  the  origin  of 
several  lines  observed  in  Si: In.  The  three  strongest  of  these  lines,  at 
1108.6,  1115.9,  and  1117.7  meV,  have  been  studied  in  detail.  In  par- 
ticular, the  temperature  dependence  and  the  time  dependence  of  these 
lines  have  been  measured. 

A.  SPECTRUM  OF  Si: In 

Figure  5 gives  a spectrum  for  a Si: In  sample  grown  by  the 
Czochralski  method.  The  spectrum  shows  not  only  the  now  well  established 
bound  exciton  (BE)  lines,  but  also  several  lines  of  unknown  origin. 

B.  DECAY  MEASUREMENTS  ON  THE  A,  B,  C LINES 

The  decay  transients  of  the  A,  B,  C lines  have  been  measured  as 

a function  of  temperature.  Typical  results  for  15  K are  shown  in 

Figure  6.  The  most  striking  aspect  of  these  decay  transients  is  the 

very  large  value  for  the  decay  time,  several  hundred  microseconds. 

These  times  should  be  compared  with  those  obtained  for  other  lines.  For 

example,  the  decay  time  for  free  excitons  in  the  purest  Si  available  is 
13 

about  2.6  Usee,  end  the  decay  time  for  excitons  bound  to  neutral 

14 

impurities  is  less  than  a few  ysec.  Hence,  the  decay  times  of  A and 
B are  approximately  two  orders  of  magnitude  larger  than  the  decay 
times  usually  observed  with  this  measurement  technique. 

These  very  long  decay  times  suggest  that  the  center  that  is  binding 
the  exciton  does  not  contain  a free  carrier.  For  if  the  center  did  con- 
tain a free  carrier,  the  Auger  decay  rate  involving  the  ionization  of  this 
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INTENSITY,  ARBITRARY  UNITS 


T = 20  K 7569  9 


Figure  5. 

The  luminescence  from  a Czochralski-grown  Si: In  sample. 
The  lines  labeled  BE  are  due  to  an  electron-hole  pair 
bound  to  a neutral  In.  The  lines  labeled  with  A,  B, 

C,  and  ? are  of  unknown  origin. 


Figure  6. 

The  decay  characteristics  of 
the  lines  labeled  A and  C in 
Figure  5. 
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center  would  be  many  orders  of  magnitude  faster  than  the  observed  decay 
rate.  Hence,  these  lines  might  be  interpreted  as  resulting  from  excitons 
bound  to  isoelectronic  traps.  Since  these  lines  have  only  been  observed 
in  samples  doped  with  In,  this  isoelectronic  trap  is  expected  to  be  a 
complex  involving  perhaps  an  In  and  some  other  species  (e.g.,  vacancy, 

0,  N,  C) . These  lines  are  the  first  published  results  on  the  decay  of 
isoelectronic  centers  in  Si. 

C.  TEMPERATURE  DEPENDENCE 

To  examine  the  question  of  whether  any  one  of  the  lines  A,  B,  or  C 
involves  different  states  of  excitons  bound  to  a center,  we  measured  the 
temperature  dependence  of  the  ratio  of  the  intensities  of  these  lines. 

If  two  of  the  lines  result  from  different  initial  states  of  the  same 
center,  then  the  ratio  should  vary  as 

-AE/kT 
e , 

where  AE  is  the  splitting  between  these  lines.  The  variation  of  these 
ratios  is  plotted  in  Figure  7.  The  fact  that  none  of  these  ratios  has 
the  expected  dependence  indicates  that  they  do  not  result  from  different 
initial  states  on  the  same  center.  The  data  show  that  the  ratio  of  A to 
B is  independent  of  temperature,  which  might  suggest  that  these  lines 
occur  because  of  different  final  states  of  the  impurity.  However,  since 
the  A and  B lines  have  different  decay  times,  they  do  not  occur  because 
of  transitions  resulting  from  the  same  initial  state  going  to  different 
final  states. 

The  results  of  measurements  of  the  temperature  dependence  of  the 
decay  times  are  given  in  Figure  8.  These  data  show  a temperature  depend- 
ence that  is  at  least  partially  understood.  If,  at  high  temperatures, 
a carrier  or  an  exciton  is  released  from  the  center,  this  would  explain 
the  observed  decrease  in  the  decay  time.  Plots  of  the  high-temperature 
data  versus  inverse  temperatures  yield  values  of  the  ionization  energy 
of  29.4  ± 0.8  meV,  29.5  ± 0.4  meV,  and  28.3  ± 0.2  meV  for  the  A,  B,  C 
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INTENSITY  RATIO 


Figure  8.  The  variation  of  the  decay  times  of  the 


lines  labeled  A,  B,  and  C in  Figure  5» 
as  a function  of  temperature. 


lines,  respectively.  These  values  are  about  10  meV  smaller  than  the 
binding  energy  of  the  exciton  obtained  by  assuming  that  the  observed 
line  is  due  to  a non-phonon  transition.  However,  it  is  possible  that 
the  ionization  is  of  a single  carrier  rather  than  of  the  electron  and 
hole  into  an  exciton  state. 

D.  SUMMARY 

The  observation  of  the  long  decay  times  for  the  A,  B,  and  C lines 
indicates  that  they  cannot  be  directly  associated  with  the  X level.  The 
long  decay  times  suggest  an  isoelectronic  center,  and  we  know  that  the 
X level  is  electrically  active  (i.e.,  acts  as  a true  acceptor).  The 
small  amount  of  information  we  have  on  the  relative  intensity  of  these 
lines  with  respect  to  the  bound  exciton  on  the  In  as  a function  of  anneal 
treatment  suggests  that  they  increase  in  number  when  the  X level  increases 
and  decrease  when  the  X level  decreases. 

The  present  activity  is  aimed  at  exploring  the  variation  of  the 
other  lines  of  unknown  origin  with  anneal  treatment.  This  program  is 
aimed  at  identifying  which  line  is  due  to  the  X level. 
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SUMMARY 

Progress  during  this  period  consisted  of  efforts  aimed  at 
characterizing  the  Si: In  X level.  Several  measurement  techniques  were 
used  to  search  for  correlations  between  the  X-level  concentration  and 
other  crystal  defects  or  growth  conditions.  These  included  IR  trans- 
mission studies  to  determine  oxygen  and  carbon  concentrations  and  growth 
studies  to  determine  the  effects  of  gas  overpressure;  melt  surface  to 
volume  ratio;  and  crucible,  polysilicon,  and  doping  materials. 

Other  studies  continuing  during  this  period  included  investigations 
of  new  characterization  techniques  to  verify  the  X-level  determinations 
made  by  the  usual  Hall-effect  and  photoconductive  techniques.  These 
new  techniques,  photoluminescence  and  MOSFET  studies,  may  lead  to  new 
insights  into  the  nature  of  the  X level.  For  example,  a large  number  of 
lines  are  observed  in  the  photoluminescence  of  Si: In.  These  lines  may 
contain  spectral  information  about  the  X level.  In  addition,  these 
techniques  may  become  useful  general  analytical  tools  for  semiconductor 
characterization. 

The  most  significant  result  during  this  period  was  the  discovery 

E 

of  a correlation  between  the  X-level  concentration  and  the  concentration 
of  carbon.  The  major  thrust  of  the  program  will  be  to  confirm  this 
result  and  to  verify  that  growth  in  a carbon-free  environment  will  yield 
low  (or  zero)  X level  Si: In. 
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ABSTRACT 

The  properties  of  indium-doped  silicon  as  a material  for  use  as  an 
extrinsic  photoconductor  was  investigated.  Detector  structures  were 
fabricated  from  characterized  material.  Responsivi ty,  D*,  current- 
voltage  characteristics,  and  spectral  photoconductive  response  were  ob- 
tained as  functions  of  bias  and  temperature.  Carrier  lifetime  and  recom- 
bination coefficient  will  also  be  presented  as  functions  of  bias,  tempera- 
ture, and  donor  concentration.  The  effect  of  the  X level  on  the  maximum 
operating  temperature  for  various  photon  backgrounds  will  be  discussed. 


1.  INTRODUCTION 

We  are  studying  Sirin  material  for  use  in  infrared  detectors  to  determine  the  suitability  of 
these  crystals  for  monolithic  focal  plane  array  (MFPA)  applications.  Discrete  detectors  are  fabri- 
cated and  studied  to  determine  material  properties  and  to  correlate  with  Hall-effect  data. 

The  material  was  grown  by  both  float-zone  and  Czochralski  methods  in  the<lll>  and  <100> 
orientations.  Two-inch-diameter  Czochralski  ingots  were  grown  in  boch  large  and  small  melt  volume 
growers.  Detectors  were  cut  from  areas  of  wafers  adjacent  to  Hall-effect  samples  so  that  accurate 
comparisons  could  be  made  between  measurements  made  on  the  two  types  of  samples.  Detector  contacts 
were  formed  by  ion  implantation  of  boron  and  by  chrome-gold  metallization.  They  were  mounted  in  a 
variable-temperature  dewar  that  had  cold  germanium  and  room-temperature  KRS-5  windows.  The  300  K 
background  was  limited  to  a 30°  field  of  view  by  a cold  aperture.  The  signal  source  for  responsivity 
measurements  was  a 500  K blackbody  chopped  at  1000  Hz.  The  signals  and  noise  were  amplified  by  a 
room  temperature  transimpedance  amplifier  with  an  adjustable  feedback  resistor  and  measured  with  a 
tuned  voltmeter.  Photoconductive  spectral  responses  were  taken  with  a modified  Beckman  IR-12 
spectrometer . 

2.  RESPONSIVITY 

Responsivity,  measured  as  a function  of  bias  voltage  and  temperature,  is  plotted  in  Figure  1 as 
a function  of  bias  current.  This  result  for  Si: In  is  the  same  as  observed  in  Si:Ga,*  that  is,  the 
responsivity  varies  linearly  with  bias  current  at  temperatures  and  fields  low  enough  that  excess 
carriers  are  not  generated,  either  because  of  thermal  ionization  or  the  onset  of  impact  ionization. 
The  departure  of  the  data  from  the  unity  slope  line  at  the  higher  temperatures  is  a direct  measure 
of  the  current  resulting  from  thermal  ionization;  the  bias  current  is  the  sum  of  this  current  and 
the  photocurrent  produced  by  the  background.  There  is  a dc  component  of  current  produced  by  the 
500  K blackbody,  but  this  contribution  is  small  compared  to  that  resulting  from  the  background. 

* 

This  research  was  supported  by  the  Defense  Advanced  Research  Projects  Agency  and  monitored  by  the 

Army  Night  Vision  Laboratory  under  Contract  No.  DAAK70-77-C-0082 . 
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Peak  responsivities  at  6.4  pm  in  excess  of  50  A/W  were  measured.  The  responsivity , normalized  to 
bias  field,  for  this  sample  and  another  is  plotted  in  Figure  2 as  a function  of  temperature  for  a 
range  of  biases  from  200  to  6500  V/ cm.  The  curves  become  progressively  flatter  (the  reason  is 
unknown)  as  the  temperature  or  bias  is  increased,  which  reduces  the  MPFA  cooler  power  requirement. 


3.  EFFECT  OF  X LEVEL  ON  OPERATING  TEMPERATURE 
The  detectivity  of  an  Si: In  extrinsic  detector  is  given  by 


D*  = 


2hc  V **  + ^TH 


(1) 


where  the  equivalent  thermal  background  flux  is  given  by 


a(n)  BN 


In 


kT 


TH 


(2) 


Eq.  2 holds  only  if  In  dominates  the  thermal  equilibrium  carrier  concentration  at  low  temperatures. 


If  the  X level  dominates  this  behavior,  then  Q is  given  by 
. TH 


a(n)  BNV  Nx 


_X 

kT 


TH 


°ln8  NIN 


(3) 


If  becomes  so  small  that  from  Eq.  3 becomes  from  Eq.  2,  then  Eq.  3 will  no  longer 

be  valid,  and  Eq.  2 should  be  used.  Also,  if  N^  < N^,  the  X level  no  longer  dominates  the  low- 
temperature  carrier  concentration,  and  a transition  is  made  to  Eq.  2,  which  becomes  valid  for 


NX  ^ V . 


By  setting  Q = Q , expressions  for  Tdt  can  be  derived  to  yield 
IH  B BLIP 


BLIP 


E_  /k 
In 


ot(n)BN, 


In 


°IngQB 


(4) 


if  Eq.  2 is  valid  and 


Vk 


BLIP 


In 


a(n)BNvNx 


°lnEQBNIn 


(5) 


if  Eq.  3 is  valid. 

To  test  the  a 

(N^/Nj^)]  as  a function  of  because  Eq.  5 can  be  rearranged  to  be 


To  test  the  agreement  of  the  data  with  the  theory,  it  is  convenient  to  plot  ln[(cx(n)  Ny/QB) 


a(n)Ny  Nx 


O-  g E /kT 
In°  . X BLIP 


QB  Nln 


(6) 
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Note  that  Q is  on  the  left  hand  side  of  the  equation,  in  contrast  to  the  arrangement  in  Ref,  2, 

This  is  to  account  for  small  variations  in  Q between  samples.  These  quantities  are  shown  plotted 

B 

in  Figure  3 along  with  the  best-fitting  line  of  the  expected  0.111  eV  slope.  The  large  error  bar 

shown  for  one  point  was  due  to  uncertainty  introduced  by  drift  in  that  sample.  The  error  bar  is  not 

typical  of  the  other  samples  for  which  significantly  smaller  errors  are  estimated.  The  hatched 

area  on  the  left  side  of  the  figure  shows  the  range  of  value  of  Tg^ip  ^ounc*  ^rom  E<1.  4 (that  is, 

assuming  that  X has  no  effect  for  the  samples  plotted  in  the  figure).  The  theoretical  curve  should 

make  an  asymptotic  transition  between  the  0.111-eV  line  and  the  vertical  hatched  area  (there  is  a 

family  of  transitions,  slightly  different  for  each  sample).  The  lowest  point  shown  in  Figure  3 

falls  somewhat  below  the  0.111-eV  line  and  probably  lies  on  such  a transition  curve.  T for 

this  sample  is  only  2.8°  below  the  maximum  possible  value  of  T calculated  from  Eq.  4 for  this 

BLlr 

sample  and  background. 

—6  3 

The  data  is  in  excellent  agreement  with  the  theoretical  line,  drawn  for  B * 1.2  x 10  cm  /s. 

*“7  “6  3 

This  value  falls  well  within  the  range  of  values  for  B of  4 x 10  to  4 x 10  cm  /s  reported  in  the 
3 

literature. 

The  temperature  dependence  of  D*  for  several  Si: In  detectors,  both  discrete  and  monolithic,  is 

plotted  in  Figure  4.  The  data  on  the  discrete  detectors  (the  solid  symbols)  were  all  taken  at  the 

same  bias  and  background  conditions.  Thus,  the  differences  in  the  D*  rolloff  temperature,  T_. , 

BLlr 

are  due  to  differences  in  the  ratio  of  Nv  to  N_  , as  discussed  above.  Sample  C64701A  was  made  from 

X In 

material  grown  in  the  large  Czochralski  furnace  at  Hughes  Industrial  Products  Division  (Carlsbad, 
California).  This  sample  had  a significantly  lower  X-level  concentration  than  crystals  grown  in  the 
smaller  HRL  furnaces.  Also  shown  in  the  figure  are  data  taken  at  several  backgrounds  from  a mono- 
lithic discrete  detector. 

Theoretical  curves  are  shown  of  D*  for  the  no  X-level  case  and  for  several  N„  to  NT  ratios. 

X In 

These  curves  were  obtained  by  solving  the  T equations  given  above  as  a function  of  Qn,  which  was 
^ BLIP  B 

converted  to  D . The  right-most  solid,  no-X  curve  gives  the  optimum  for  Si: In  photoconductive 
BLlr 

detectors.  The  dashed  curves  show  the  degradation  resulting  from  various  amounts  of  X.  The  dashed 

curves  follow  a 0.111-eV  slope,  and  the  solid  curve  follows  0.1572-eV  slope.  As  the  background  is 

decreased,  a given  N /N  imposes  a relatively  more  severe  penalty.  Thus,  while  sample  C64701A 
x in 

approaches  pure  In  performance  at  Q = 3 x 10  , the  small  amount  of  X in  the  sample  may  still 

B 11 

impose  a 2°  or  3°  temperature  penalty  at  Q = 7 x 10  . Of  course,  this  represents  a significant 

B 

improvement  over  the  previous  high-X  samples. 


4.  CARRIER  LIFETIME 

The  carrier  lifetime  as  a function  of  temperature  and  bias  field  was  calculated  from  respons- 
ivitv  measurements.  The  usual  equation  for  responsivity  as  a function  of  quantum  efficiency  and 
photoconductive  gain  can  be  rearranged  to  solve  for  lifetime: 

T . Rhv  i!  > (7) 

en  pv 

where  R is  the  measured  responsivi ty , hv  is  the  photon  energy,  e is  the  electronic  charge,  n is  the 
quantum  efficiency  (which  can  be  obtained  either  from  D*  measurements  or  from  measurements  of  In  con- 
centration and  absorption  coefficient),  Z is  the  sample  thickness  (electrical  length),  p is  the 

2 

R.  Baron,  M.H.  Young,  J.F.  Baukus,  O.J.  Marsh,  and  J.J.  Sheets,  Proc.  IRIS  Det.  Spec.  Group  Meet. 
(1977). 

*A.C.  MiLnes,  Deep  Impurities  in  Semiconductors  (John  Wiley  and  Sons,  New  York,  1973),  p.  290. 
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carrier  mobility,  and  V is  the  bias  voltage.  In  this  case,  n is  calculated  from: 


n 


o N l 
A In 


(8) 


where  a A is  the  absorption  cross  section,  and  is  the  indium  concentration  measured  by  the  Hall 

effect.  The  mobility  used  in  Eq.  7 is  calculated  to  include  both  the  neutral  scattering  effects  of 

A 5 

the  indium  doping  and  the  high-yield  velocity  saturation  effects  as  measured  by  Ottaviani  et  al. 


U 


(9) 


where  y is  the  value  obtained  from  low-field  Hall-effect  measurements,  and  y , the  hot  carrier 

HALL  HUi 

mobility  from  Ref.  4,  is  a function  of  bias  field,  temperature,  crystal  orientation,  and  conductivity 
type. 

The  carrier  lifetime  is  a function  of  the  compensating  dopant  concentration  and  the  recombina- 
tion coefficient,  B: 


(10) 


Figure  5 is  a plot  of  T as  obtained  from  responsivity  measurements  versus  (N  - N ) as  measured  by 
the  Hall  effect  for  several  bias  fields.  The  lines  represent  the  least-squares  fit  of  Eq.  10  to  the 
data,  and  the  resulting  B values  are  indicated.  The  data  represent  samples  grown  by  the  float-zone 
and  Czochralski  techniques  in  both  large  and  small  furnaces.  Also,  the  x and  y values  were  obtained 
by  completely  independent  measurements. 

As  the  bias  field  is  increased,  the  carrier  lifetime  increases.  This  is  consistent  with  the 
Lax  giant  trap  model, ^ which  states  that  carrier  recombination  is  through  high  excited  states  of  the 
dopant.  As  the  electric  field  is  increased,  lowering  the  barrier,  holes  in  the  excited  states  have 
an  increasing  probability  of  returning  to  the  valence  band  and  thus  "live"  longer  as  carriers.  This 
barrier  lowering  is  the  Poole-Frenkle  effect. 

These  data  were  calculated  from  measurements  taken  at  40,  45,  and  50  K for  fields  of  approxi- 
mately 200,  400,  and  1800  V/cm.  Fields  varied  from  sample  to  sample  by  about  10%,  which  was  not 
fully  accounted  for  in  the  calculations.  The  values  of  B obtained  by  this  procedure  for  these 
temperatures  and  fields  are  plotted  in  Figure  6 as  a function  of  temperature.  As  temperature  is 
increased,  B becomes  smaller  (T  increases).  This  is  because  the  larger  values  of  kT  encompass  more 
and  more  excited  states,  which  increases  the  probability  for  captured  holes  to  return  to  the  band. 

g 

This  effect  is  more  fully  discussed  in  another  paper.  B also  decreases  with  increasing  electric 
field,  as  discussed  above.  The  data  show  good  correlation  with  both  temperature  and  field,  although 
more  data  over  a wider  range  of  the  parameters  are  needed. 
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"*G.  Ottaviani,  L.  Reggiani,  C.  Canali,  F.  Nava,  and  A.  Alberigi-Quaranta,  Phys.  Rev.  B12,  3318  (1975). 
^M.  Lax,  J.  Phys.  Chem.  Solids  8,  66  (1959). 

7J.L.  Hartke,  J.  Appl.  Phsy.  39,  4871  (1968). 
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The  hole  capture  cross  section,  o,  can  be  calculated  from  the  recombination  coefficient: 

B “ vnf  • 

where  B is  the  hole  capture  coefficient,  and  is  the  thermal  equilibrium  carrier  velocity: 


TH 


/8kr\I/2 


L TTm 


n 


(id 


(12) 


where  m*  is  the  hole  effective  mass.  Although  some  form  of  the  conductivity  effective  mass  would  be 

9 

more  appropriate,  the  density  of  states  effective  mass  is  used  because  it  is  known.  The  resulting 
discrepancy  is  not  large. 

The  calculated  capture  cross  section  is  plotted  in  Figure  7 as  a function  of  temperature.  Also 
plotted  are  the  results  of  Norton  et  al.^°  on  Si:P  at  lower  temperatures  and  a range  of  values  for 
Si:ln  at  77  K reported  by  Milnes. * Our  data,  which  appear  to  be  in  quite  good  agreement  with  these 
values,  provide  information  on  the  field  dependence  of  the  cross  section.  The  agreement  with  the 
Si:P  data  is  to  be  expected  since  the  higher  excited  state  spectra  are  similar  for  all  dopants  in 
silicon. 

5.  SUMMARY 

We  have  measured  the  temperature  and  bias  dependence  of  Si: In  material  from  ingots  that  had 
been  used  to  produce  successful  MFPAs.  We  found  that  the  responsivity  is  directly  proportional  to 
the  bias  current  for  the  usual  fields  and  temperatures  of  MFPA  operation.  The  responsivity 
becomes  less  temperature  sensitive  at  higher  fields  and  temperatures.  The  results  of  a refined  cal- 
culation of  the  effect  of  the  X level  on  detector  performance  were  compared  with  experiment.  The 
amount  of  X level  in  the  present  material  causes  a temperature  penalty  of  only  2°  to  3°  at  low  back- 
grounds. We  have  also  calculated  the  field  and  temperature  dependence  of  the  carrier  recombination 
coefficient  and  cross  section.  To  our  knowledge,  this  is  the  first  time  that  field-dependence  data 
have  been  presented. 
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FIGURE  U.  D*  VERSUS  TEMPERATURE.  Peak  D* 
versus  temperature  for  Si: In  discrete  and  mono- 
lithic discrete  detectors.  Also  shown  are  theo- 
retical curves  for  several  to  N ratios. 
Photon  background  units  are  photons?sec-cm^. 


FIGURE  3.  EFFECT  OF  X LEVEL  ON  TfiL  . The  effect 
of  the  Nj^  to  NIn  r*.tio  on  the  D*  rolloff  temper- 
ature, TBLip>f°r  several  samples.  The  curve  is  a 
fit  to  the  data  using  Eq.  6 with  B as  the  fitting 
parameter.  The  shaded  area  on  the  left  is  the 
range  of  TgLIP  values  that  would  be  obtained  if 
the  samples  had  no  X level. 
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FICURE  7.  CAPTURE  CROSS  SECTION  VERSUS 
TEMPERATURE.  Hole  capture  cross  section 
as  a function  of  temperature  for  Si: In 
at  three  bias  fields.  The  range  of 
values  indicated  at  77  K are  from  Ref.  3. 
The  values  at  10  to  20  K.  are  for  Si:P 
and  are  from  Ref.  10. 
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PHOTOLUMINESCENCE  STUDIES  OF  ACCEPTORS  IN  SILICON* 

S.A.  Lyon,  K.R.  Elliott,  G.  Mitchard,  D.L,  Smith,  and  T.C.  McGill 

California  Institute  of  Technology 
Pasadena,  California  91125 

and 

J.P.  Baukus,  R.  Baron,  M.H,  Young,  and  O.J.  Marsh 

Hughes  Research  Laboratories 
Malibu,  California  90265 

ABSTRACT 

The  photoluminescence  spectra  of  doped  Si  are  summarized,  with 
particular  emphasis  on  the  spectra  caused  by  the  acceptors.  The  processes 
governing  the  relative  intensity  of  the  luminescence  lines  are  discussed, 
and  data  is  presented  that  makes  it  possible  to  estimate  the  relative 
intensity  of  the  lines.  New  results  are  reported  for  the  luminescence 
spectra  of  Si: In,  including  the  observation  of  some  unusually  long-lived 
lines. 

1.  INTRODUCTION 

Photoluminescence  from  excitons  bound  to  impurity  centers  (BE)  has  been  used  as  a diagnostic 
technique  for  several  years.  The  primary  emphasis  of  the  work  has  been  the  III-V  and  II-VI  semicon- 
ductors, which  have  potential  as  optoelectronic  devices.^  In  the  last  few  years  a great  deal  of 

interest  has  developed  in  the  luminescence  spectra  of  silicon,  both  because  it  is  an  ideal  system  in 

2-5 

which  to  study  the  fundamentals  of  the  luminescence  process  and  because  of  potential  for  using 
the  luminescence  spectrum  as  a diagnostic  measure  of  the  impurity  and  defect  concentration.  This 
paper  summarizes  some  of  our  recent  results  on  the  luminescence  spectra  of  the  acceptors  in  silicon 
with  particular  emphasis  on  Si: In. 

Section  2 briefly  describes  the  experimental  technique.  Section  3 presents  some  of  the  exper- 
imental data  on  the  line  positions  for  B,  P,  Al,  Ga,  In,  and  Tl,  plus  some  of  the  complications  of 
the  spectra.  Section  4 presents  the  data  for  two  of  the  most  important  processes  governing  the 
photoluminescence  process:  the  capture  of  an  exciton  on  the  impurity  and  the  predominant  decay 
mechanisms  of  the  exciton  once  it  is  bound  to  the  impurity.  Section  5 presents  the  preliminary 

This  research  was  supported  by  the  Defense  Advanced  Research  Projects  Agency  and  monitored  by  the 
Army  Night  Vision  Laboratory,  under  Contract  No.  DAAK70-77-C-0082. 

1See,  for  example,  A. A.  Bergh  and  P.J,  Dean,  Light  Emitting  Diodes  (Clarendeon  Press,  Oxford,  1976), 
and  the  references  contained  therein, 

2S.A.  Lyon,  D.L.  Smith,  and  T.C.  McGill,  Phys.  Rev.  B17,  2620  (1978). 

3K.R.  Elliott,  G.C.  Osbourn,  D.L.  Smith,  and  T.C.  McGill,  Phys.  Rev.  B17,  1808  (1978). 

4S.A.  Lyon,  G.C.  Osbourn,  D.L.  Smith,  and  T.C.  McGill,  Solid  State  Comm.  23.,  425  (1977). 

5K.R.  Elliott,  D.L.  Smith,  and  T.C.  McGill,  Solid  State  Comm.  24,  461  (1978). 
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results  from  a study  of  the  luminescence  from  Si: In  samples;  we  have  observed  several  previously 
unknown  lines,  some  with  unusually  long  lifetimes.  Section  4 summarizes  our  conclusions. 

2.  EXPERIMENTAL  PROCEDURES 

The  experimental  setup  is  shown  in  Figure  1.  The  Si  samples  were  prepared  by  lapping  and  then 
etching  with  HN0^:HF(7:1)  and  cleaning  with  methanol.  The  samples  were  then  placed  in  a helium  cryo- 
stat, and  the  bath  temperature  was  measured  using  a Ge  temperature  sensor  in  the  sample  block. 
Electron-hole  pairs  were  produced  by  irradiating  the  sample  with  above  bandgap  radiation  from  a GaAs 
or  argon  ion  laser.  The  luminescence  spectrum  resulting  from  the  recombination  of  an  electron-hole 
pair  bound  to  a neutral  impurity  was  analyzed  with  a spectrometer,  detected  by  a cooled  S-l  photo- 
multiplier, and  processed  using  standard  photon  counting  techniques.  The  type  of  signal  observed  is 
illustrated  in  the  figure.  It  consists  of  a series  of  peaks  of  rather  widely  varying  Intensities. 

3.  LUMINESCENCE  SPECTRA  OF  THE  DOPANTS  IN  SILICON 

The  luminescence  spectrum  of  a sample  of  Si:Al  is  shown  in  Figure  2.  This  data  shows  the  spec- 
trum in  the  region  where  the  recombination  occurs  without  a phonon,  the  no-phonon  region.  Other 
replicas  of  this  spectrum  are  observed  at  energies  corresponding  to  this  spectrum  except  shifted  by 
a phonon  energy,  the  so  called  phonon-assisted  transitions.  The  spectrum  consists  of  several  lines, 
some  of  them  due  to  the  A1  dopant  and  others  due  to  the  background  P impurity  in  the  sample.  All  the 
lines  labeled  P are  due  to  the  impurity  phosphorus.  The  line  at  approximately  1.1499  eV  is  due  to 
the  decay  of  a single  electron-hole  pair  bound  to  P.^  The  other  lines  labeled  P are  due  to  the 
decay  of  a single  electron-hole  pair  bound  to  a center  with  varying  numbers  of  other  electron-hole 
pairs  bound  to  it.  The  lines  labeled  BE  are  due  to  the  decay  of  a single  electron-hole  pair  bound 
to  a neutral  Al.  The  other  lines,  labeled  b^,  b and  b^,  are  due  to  the  decay  of  an  electron-hole 
pair  bound  to  a.  neutral  center  along  with  several  other  electron-hole  pairs. ^ 

This  spectrum  illustrates  three  important  features.  First,  both  donors  and  acceptors  can  be 
observed  simultaneously  in  luminescence.  Second,  the  lines  for  different  impurities  are  displaced 
from  each  other  in  such  a way  that  they  can  be  resolved  spectrally.  Third,  each  impurity  produces 
several  lines,  which  give  a unique  signature  for  this  impurity. 

Figure  3 presents  positions  of  the  strongest  line  for  each  impurity,  the  so  called  principal 
bound  exciton  line.  Again,  this  data  illustrates  the  uniqueness  of  the  spectral  position  of  the 
line.  Further,  the  energy  positions  of  the  lines  for  a given  impurity  are  ordered  according  to  the 
strength  of  the  binding  of  the  neutral  impurity.  This  observation  has  been  quantified  to  give  the 
approximate  Haynes  rule.^ 


4.  THE  PHOTOLUMINESCENCE  PROCESSES 

To  understand  and  quantify  the  observed  photoluminescence  spectrum  requires  an  understanding  of 
the  various  parts  of  the  photoluminescence  processes.  To  obtain  information  on  the  relative  contri- 
butions of  the  various  impurities,  we  need  to  understand  the  capture  of  the  free  exciton  on  the 


^K.  Kosai  and  M.  Gershenzon,  Phys.  Rev.  B9,  723  (1974). 
^J.R.  Haynes,  Phys.  Rev.  Lett.  4_,  361  (1960). 
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neutral  impurity  and  the  processes  governing  the  relative  strength  of  the  various  lines.  The 
capture  process  is  one  in  which  a free  exciton  is  captured  by  a neutral  impurity,  as  characterized 
by  a capture  cross  section.  The  measured  capture  cross  section  for  Si: In  is  presented  in  Figure  4. 

The  data  shows  that  the  capture  cross  section  varies  very  rapidly  with  temperature:  it  is  about 
-15  -13  2 

10  cm  at  30  K and  ibout  10  cm  at  10  K.  The  data  shows  a variation  somewhat  similar  to  that 
found  for  the  capture  of  a carrier  on  an  ionized  center,  where  the  mechanism  is  capture  into 

g 

highly  excited  states  of  the  impurity  center,  the  Lax  giant  trap  mechanism.  Since  we  expect  the 

highly  excited  states  of  the  bound  excitons  to  be  very  similar  for  different  impurities,  this  cross 
section  data  should  also  apply  to  other  impurities. 

To  estimate  the  relative  strength  of  lines  resulting  from  various  impurities,  we  assume  that 
the  capture  cross  sections  are  the  same  for  the  impurities  and  that  pump  intensity  is  such  that  the 
lines  are  not  saturated.  That  is,  we  assume  that  most  of  the  impurity  centers  are  not  occupied  by 
electron-hole  pairs.  With  these  reasonable  assumptions,  the  relevant  factors  in  the  intensity,  I,  of 
a given  line  are  given  by 


I « N 


"nr  + *r 


(1) 


where  N is  the  number  of  centers,  and  R and  R are  the  radiative  and  nonradiative  rates,  respect- 

R NK 

ively.  For  Si,  we  can  write 


This  can  be  used  to  rewrite  Eq.  1 as 


I a N 


(2) 


The  radiative  rate  is  proportional  to  the  oscillator  strength,  f,  for  a given  transition: 


f . (3) 

In  Figure  5,  we  have  plotted  the  oscillator  strength  for  the  acceptors  in  the  no-phonon  transition. 
The  data  shows  that  oscillator  strength  increases  with  the  binding  of  the  acceptors  according  to 
the  approximate  relation 


f “ . (4) 

where  E is  the  binding  energy  of  the  hole  on  the  acceptor.  The  predominant  decay  process  for 

A 9 

excitons  on  these  impurities  in  Si  is  the  nonradiative  Auger  process.  In  this  process,  the 

electron-hole  pair  on  the  center  recombines  and  the  energy  is  carried  away  by  the  remaining  carrier 

on  the  impurity  center.  The  lifetime  of  the  bound  excitons  on  B,  Al,  and  Ga  due  to  this  process  is 

shown  in  Figure  6.  Schmid^  has  made  measurements  of  the  lifetime  for  In,  P,  and  A.  He  found 

V Lax,  Phys.  Rev.  119,  1502  (1960). 

^G.C.  Osbourn  and  D.L.  Smith,  Phys.  Rev.  B16,  5426  (1977). 


47 


UNCLASSIFIED 


UNCLASSIFIED 


T_  = 2.7  nsec,  T * 272  nsec,  and  T ■ 183  nsec.  The  Auger  decay  rate  Increases  rapidly  with 

An  v 9 io  As  9 10 

increasing  binding  energy.  * This  increase  is  given  approximately  by  ' 


(5) 


where  n is  between  4 and  6. 

Combining  the  results  of  Figures  5 and  6 or  the  approximate  dependences  given  by  Equations  4 and 
5,  we  can  obtain  the  correction  relating  the  relative  intensities  to  the  relative  number  of  cen- 
ters. Using  the  approximate  expression,  the  correction  factor  varies  as  E m,  where  m is  in  the 
range  of  1 to  3.  Thus,  we  find  that  the  intensity  for  a given  concentration  decreases  with  increas- 
ing binding  energy  of  the  impurity  but  not  as  fast  as  the  Auger  rate  would  suggest. 


5.  PHOTOLUMINESCENCE  SPECTRA  OF  Si: In  SAMPLES 

The  interest  in  Si: In  for  making  infrared  detectors  and  the  need  to  understand  the  origin  of 
the  "x-level"^  has  led  to  a study  of  the  luminescence  of  Sirin  samples.  In  Figure  7 is  the  lumin- 
escence spectrum  of  one  of  the  Sirin  samples.  The  spectrum  is  dominated  by  several  very  intense 
lines  (labeled  BE)  that  are  due  to  the  In  dopant.  Other  lines  (labeled  A,  B,  and  C)  are  intense  and 
of  unknown  origin.  They  have  not  been  observed  in  samples  doped  with  Ga  or  A1  or  any  of  the  donors. 
These  lines  tend  to  be  much  stronger  in  Czochralski-grown  material  as  compared  to  samples  prepared 
by  the  float  zone  technique.  At  one  point,  these  lines  were  thought  to  be  likely  candidates  for  the 
Mx-level”  (see  the  discussion  of  the  transients  below).  Several  other  lines  of  unknown  origin  are 
also  indicated  in  the  spectrum  (by  question  marks).  The  positions  of  all  of  the  lines  and  their 
assignments  are,  where  possible,  given  in  Table  1. 

The  lines  labeled  A,  B,  and  C were  investigated  by  studying  their  sample  dependence,  decay 
characteristic,  and  temperature  dependence.  The  most  unusual  property  of  each  line  is  its  decay 
characteristic.  Figure  8 shows  the  decay  characteristics  of  the  A and  C lines.  The  measured  decays 
are  exponential  and  yield  decay  times  on  the  order  of  a few  hundred  microseconds.  These  decay 
times  are  unusually  long  for  excitons  in  Si.  All  of  the  other  decay  times  measured  for  bound  exci- 
tons  in  Si  are  typically  less  than  a few  microseconds.  Hence,  these  lines  possess  decay  times  that 

are  at  least  two  orders  of  magnitude  larger  than  other  lines  observed  in  Si.  Very  recently,  Weber, 
12 

Schmid,  and  Sauer  observed  lines  with  different  energies  than  those  reported  here  but  with  very 

long  decay  times.  They  interpret  these  lines  as  resulting  from  isoelectronic  traps  associated  with 
carbon  in  Si.  It  seems  very  likely  that  the  lines  we  have  observed  are  due  to  an  isolectronic  cen- 
ter in  the  Si  and  that  the  observed  decay  time  is  the  radiative  decay  time  for  the  center.  Since 
these  lines  have  only  been  observed  in  samples  doped  with  In,  it  seems  likely  that  they  involve  In. 

We  are  presently  studying  the  other  lines  in  the  Si: In  spectrum  and  hope  to  report  on  them 
shortly. 


10w.  Schmid,  Phys.  Stat.  Solid!  b84,  529,  (1977). 

liR.  Baron,  M.H.  Young,  J.K.  Neeland,  and  O.J.  Marsh,  Appl.  Phys.  Lett.  30,  594  (1977), 
12 

J.  Weber,  W.  Schmid,  and  R.  Sauer,  private  communication. 
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6.  SUMMARY 

We  have  given  a brief  summary  of  some  of  our  results  on  the  photoluminescence  spectra  of 
acceptors  in  Si.  The  results  show  that  the  photoluminescence  lines  for  an  impurity  give  a unique 
signature  for  that  impurity  and  that  we  should  be  able  to  use  this  technique  as  a diagnostic  tool  for 
looking  at  the  impurities  in  infrared  detector  materials.  We  presented  the  results  for  the  param- 
eters that  would  allow  ui  to  convert  relative  line  intensity  data  to  relative  concentrations. 

Finally,  we  presented  some  new  results  on  the  Sirin  luminescence  that  suggest  that  In  forms  a 
center  in  Si  that  acts  as  an  isoelectronic  trap.  Several  new  lines  were  reported  for  the  Sirin 
samples. 

Table  1.  The  energy  and  assignment  of  some  of  the  lines  observed 
in  the  photoluminescence  spectrum  shown  in  Figure  7 for  a Sirin  sample 


Peak 

Energy,3 

MeV 

Identification 

1077.6 

? 

1082.0 

TO  + LO  phonon  In  BE  (J  = 0) 

1085.1 

TO  + LO  phonon  In  BE  (J  = 2) 

1093.3 

? 

1098.4 

TO  + L0  phonon  FE 

1101.1 

9 

1105.7 

? 

1108.6 

A 

1110.7 

9 

1114.4 

? 

1115.9 

B 

1117.7 

C 

1121.5 

TA  phonon  In  BE  (J  = 0) 

1125.0 

TA  phonon  In  BE  (J  = 2) 

1129.2 

? 

1136.5 

? 

1140.5 

In  BE  (J  * 0) 

1143.7 

In  BE  (J  = 2) 

1154.3 

9 

aEnergies  accurate  to  about  ±0.2  meV. 
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Figure  1.  Figure  2. 

Schematic  of  the  experimental  setup  and  The  photoluminescence  spectra  of  A1  doped  Si  in  the 

the  photoluminescence  process.  energy  range  for  no-phonon  assisted  transitions.  The 

lines  labeled  BE,  b^,  b2*  and  b^  are  associated  with 
A1  impurities.  The  lines  labeled  P are  associated 
with  phosphorous  impurities  in  the  Si. 


1.105  1.110  1.140  1.145  1.150  1.155 

ENERGY.  eV 


Figure  3. 

The  position  of  the  principal  bound  exciton  line  for  various 
Impurities  in  the  no-phonon  assisted  transition.  The  positions 
for  B and  P were  taken  from  Ref.  6. 
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Figure  4. 

The  measured  temperature  dependence  of  the 
capture  cross  section  of  free  excitons  on 
neutral  indium  impurities  in  silicon.  The 
data  is  given  for  different  doping  levels 
ranging  from  1.8  x 10^  to  5.5  x 10^  cm“^. 


Figure  5. 

The  total  absolute  oscillator  strength  for 
the  lowest  known  transitions  of  an  electron- 
hole  pair  bound  to  an  acceptor.  The  data 
indicated  by  ■ were  taken  from  P.J.  Dean, 
W.F.  Flood,  and  G.  Kaminsky  (Phys.  Rev.  163, 
721  (1967)) 
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Figure  6. 

Bound  exciton  lifetimes  versus  impurity  bind- 
ing energy  for  the  four  common  acceptors  in  Si. 
The  solid  circles  are  measured  values  and  the 
squares  are  calculations  of  the  Auger  lifetime. 
The  measured  value  for  Si: In  is  from  Ref.  10. 
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Figure  7. 

The  luminescence  from  a Czochralski-grown  Si: In  sample.  The 
lines  labeled  BE  are  due  to  an  electron-hole  pair  bound  to  a 
neutral  In,  The  lines  labeled  with  A,  B,  C,  and  the  lines  label- 
ed with  ? are  of  unknown  origin. 


Figure  8. 

The  decay  characteristics  of  the  lines  labeled 
A and  C in  Figure  7. 
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